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We investigate (theoretically and numerically) the dynamics of a soliton moving in an asym- 
metrical potential well with a finite barrier. For large values of the width of the well, the width of 
the barrier and/or the height of the barrier, the soliton behaves classically. On the other hand, we 
obtain the conditions for the existence of soliton tunneling with sub-barrier kinetic energies. We 
apply these results to the study of soliton propagation in disordered systems. 
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The escaping process of a particle from a potential well, 
as that shown in Fig. 1 (barrier crossing), is a problem 
£f) ■ of great importance in almost all areas of physics |jj . 

In the case of a classical particle, this escape should 
occur over the barrier with the help of external pertur- 
bations (e.g., thermally activated barrier crossing) |Q. 
On the other hand, a quantum particle can perform tun- 
neling with certain probability p < 1 . 

In the present letter we address the question: what hap- 
pens if, in the potential well, instead of a point-like parti- 
cle we have a soliton? This is very relevant to Skyrmion 
+1j . models for nucleon physics, the motion of fluxons in long 
Josephson junctions with impurities, the dynamics of do- 
O !■ main walls in ferroelectric materials in the presence of 
inhomogeneous electric fields, and many other physical 
. ,_h , systems where the solitons move in a potential created 
by inhomogeneities and external forces 

It is well-known that a soliton can behave as a classical 
particle in some physical systems |^||. However, recently 
there has been a great interest in non-classical behaviors 
of the soliton [Q-Q. In particular, we are interested in 
extremely surprising phenomena that can occur when the 
soliton behaves as an extended object [p|^0|. 

Among these phenomena is the soliton tunneling sug- 
gested by Kalbermann in a beautiful paper ]i~5[ |. How- 
ever, this was a numerical work and the tunneling with 
sub-barrier kinetic energies was not observed in the spe- 
cific situations in which the numerical experiments were 
performed. 

In the present letter we show (theoretically and numer- 
ically) that the tunneling with sub-barrier kinetic ener- 
gies is, indeed, possible! 

As an example we consider the perturbed 4> -equation: 
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<f) = F{x). (1) 
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The external force F(x) is such that a point-like soliton 
would feel an effective potential like that shown in Fig. 
1. 

When the soliton is treated as a point-like particle, the 
zeroes of F(x) are equilibrium points fl. The zeroes xq 

(F(xq) = 0) for which > are stable equilib- 

-I X — Xq 

rium positions. In the opposite case, they are unstable. 

For our theoretical calculations we will use the force 
F{x) defined in the following way: 

F{x) = Fi(x), for x < x\ (2) 



F{x) = c, for x > x*, (3) 

where Fi(x) = \A{A 2 - l)tanh(Ba;) + \A{AB 2 - 
x* [x* > 0) is the point where Fx{x) has 

a local minimum (^ dF ^ - = 0^ , and c = Fi(x*). 

The condition |F(-oo)| = \A\A 2 - l| < ^= should 
hold for the stability of the soliton as a whole. This force 
allows us to solve the problem of soliton dynamics in a 
neighborhood of the equilibrium points ^ [l0| . For in- 
stance, the stability problem (</>(a;,t) = <j>k{x) + f(x)e xt ) 
of the equilibrium point xq = is reduced to the 
eigenvalue problem Lf = Tf, where L = —d 2 + 

( |^4 2 - i - T5 — TTT5-T ) and T = -A 2 . The eigenvalues 
of the discrete spectrum are given by r„ = — i + B 2 (A + 
2An - n 2 ) where A(A + 1) = fg^. 

Our analysis reveals that if A 2 > 1 and 4B 2 < 1, the 
force given by Eqs. (||)-(||) possesses the desired proper- 
ties, i.e. there is a zero that would correspond to a stable 
equilibrium position in a point x — —d (d > 0) and a zero 
in the point x = that would correspond to an unsta- 
ble equilibrium position and would serve as a potential 
barrier. For x > the potential is a monotonically de- 
creasing function. 

In fact, if 2i3 2 (3A 2 — l) < 1, then the soliton behaves 
classically. In this case, the soliton feels the barrier in 
the point x = 0. If the soliton is situated in a vicinity of 
point x — with zero initial velocity and with the center 
of mass in a point x < 0, it will not move to the right of 
point x — 0. 

On the other hand, if 2B 2 (ZA 2 - l) > 1, the soliton 
will move to the right, crossing the barrier even if its 
center of mass is placed in the minimum of the potential 
and its initial velocity is zero (see Fig. 2). In this case 
the soliton performs tunneling with sub-barrier kinetic 
energy! 

We should remark that this phenomenon is possi- 
ble only when the distance d between the minimum 
of the potential well and the maximum of the poten- 
tial barrier holds the inequality d < 2.17, where d = 

(-g) arccosh ( y . This can be interpreted in the 

sense that the wavelength of the soliton should be com- 
parable with the width of the potential well. 
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In this context, the work H] addresses the differences 
and similarities of soliton phenomena with those of point 
particles. Although the studied phenomena are very dif- 
ferent, we can say that the main conclusion is in con- 
sistency with our previous papers [js|-flOl and the present 
one: the soliton behaves as a particle only when the width 
of the potential makes the soliton appears as point-like. 
Otherwise, the soliton can have wave-like extended char- 
acter. 

The condition for soliton tunneling can be written in 
more physical terms. In fact, the force F(x) can be de- 
fined by the value Fo = jlim^^-oo F(x) | , the width of the 
barrier S and the local maximum F m of the force between 
the points x — —d and x = (see inset in Fig. 1). In 
these terms, the approximate condition for the existence 
of soliton tunneling is Fo > Fm e s . This inequality shows 
that greater values of Fo support the tunneling, while 
greater values of F TO and S can thwart the tunneling. 

We have performed numerical experiments with the 
force F(x) defined as in Eqs. and with many 

other functions which produce an effective potential as 
that shown in Fig. 1. We have been able to control the 
values of F , F m and S. Figure 3 shows the numerical 
experiments. The dots on the curve separates two zones: 
the zone in which the soliton tunneling is possible (upper 
zone) and the zone in which the soliton tunneling is im- 
possible (lower zone). Note that the relation Fo = Fm 6 s 
is approximately satisfied. With other forces, the results 
are qualitatively equivalent. 

Sometimes, it is convenient to see the condition for the 
existence of soliton tunneling in terms of a parameter 
that defines the potential. An important characteristic 
of the potential is the height of the potential barrier: V m . 
Figure 4 shows the relationship between F and V m while 
the bifurcation condition 2B 2 (3A 2 — l) = 1 holds. For 
points above the curve, the tunneling occurs. For points 
under the curve, the soliton behaves classically. That is, 
as was expected, the height of the potential barrier is an 
opposing factor with respect to the tunneling. 

Note that the values of the potential V(x) (or the force 
F(x)) for x <C — d can influence the tunneling. This 
is in contrast with the behavior of a point-like classical 
particle. 

We should remark directly that, in this phenomenon, 
there is no violation of the energy conservation law. In 
fact, the energy of the whole system 

-£[K3)' + K59>* 

-F{x)4> + c] dx (4) 

is conserved. In this case, the soliton does not behave 
as a classical point-like particle. The soliton, as an 
extended object, possesses a wave-mechanical behavior. 
Even when the center of mass is situated in a point of 
minimal potential energy and with zero kinetic energy, 
the system as a whole can have enough energy to make 
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the soliton able to cross the barrier (although the center 
of mass goes through the barrier). Kalbermann's exam- 
ple of a high jumper is felicitous pa . However, we should 
notice the presence of nonlocal effects during this process. 

The propagation of solitons in disordered media has 
been studied intensively in last years |lq| . There is con- 
sensus in the conclusion that nonlinearity can modify the 
effects of localization and the transmission is improved. 
Nevertheless, even in a nonlinear system supporting soli- 
tons, if the latter behave as point-like particles, they can 
be trapped in the zeroes of the force F(x). We stress 
that the phenomenon of soliton tunneling can enhance 
even more the transmission. 

Consider Eq. (]l|) with F(x) defined in such a way that 
it possesses many zeroes, maxima and minima (see Fig. 
5). This system describes an array of inhomogeneities. 
The array can be studied as a series of elements with two 
zeroes and a maximum. If for each element the condi- 
tion Fo > F g is satisfied, then the soliton can cross 
the whole inhomogeneous zone (we have checked this nu- 
merically). The array can be completely disordered. If 
the condition is fulfilled, there is no localization. 

In many systems |l6|-pc| the solitons play the role of 
means of transport: they can carry energy and/or charge. 
Our result shows that, when soliton tunneling is possible, 
the soliton can be a very efficient carrier. 

Kalbermann |D| investigated impurities that are intro- 
duced in the Hamiltonian density in the following way: 

where p = p + U(x); po is a constant, U(x) is the per- 
turbation that describes the impurity. 

The perturbations used in the numerical experiments 
in Ref. are given by the function 

hi ho 

U ^ = 77 \ + 77 V ^ 

cosh 2 f 2^fiJ cosh 2 f 

This perturbation is equivalent to an effective potential 
with a maximum, a minimum or a combination of both. 
Out of the inhomogeneous zone, the potential tends to 
zero exponentially. Suppose that the soliton moves from 
the left with a kinetic energy less than the maximum of 
the potential barrier Jl5). In this case the soliton tun- 
neling does not exist. As we have shown, for the soliton 
tunneling (among other conditions) it is necessary to have 
a soliton moving in a potential well where V(x) takes val- 
ues (out of the potential well) greater than that of the 
barrier. Nevertheless, with localized impurities like the 
ones introduced in Eq. (||) , the soliton tunneling can also 
be observed. This is possible with a perturbation U (x) 
with the features shown in Fig. 6. Even in this case, 
if the soliton moves from the left (in zone A) with sub- 
barrier kinetic energy, then the soliton tunneling does not 
occur. For the tunneling, the soliton should be placed in 
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zone B. Of course, the rest of the conditions should be 
satisfied. 

The sine-Gordon soliton usually is thought to be a 
very point-like object. In fact, the unperturbed sine- 
Gordon equation is integrable and its soliton solution 
does not have discrete, internal (shape) modes. Never- 
theless, when the perturbations are not given by Dirac's 
(^-functions, the sine-Gordon soliton also is able to excite 
a great number of shape modes and behaves as an 
extended object. 

Consider the perturbed sine-Gordon equation 

4>tt - 4>xx + sin 4> = F(x), (7) 

where 

! — F for x < xq, 

_a__ Fo * for -xq < x < x , (8) 
— Fq for x > xq; 

xq and Fq are chosen such that the function F(x) is 
continuous. 

Let us see only two examples. Let Fo = 0.25 and b — 
0.55 (fixed). For a = 0.45, the soliton remains trapped in 
the potential well. It behaves classically. For a = 0.35, 
the soliton escapes from the potential well, crossing the 
barrier. We should emphasize that in both cases the force 
F(x) would correspond to a system with a potential well 
and a barrier if the soliton behaves classically. 

We would like to remark that, in our study, the posi- 
tioning of the soliton at a certain point is done consid- 
ering the equation of motion. The initial configuration 
is always a solution of the equation of motion (we can 
use both the static and the time dependent solutions in 
dependence on the physical process that led to the given 
situation). That is, we never use an initial configuration 
which the soliton can not reach by any means. For ex- 
ample, let us explain a physical situation in which the 
soliton can be positioned inside the potential well de- 
picted in Fig. 1. Suppose a soliton is captured by a 
localized inhomogeneity (there are many such situations 
discussed in the review paper || and the experimental 
papers quoted therein). Then, we apply an external con- 
stant force (e.g. for the Joscphson fluxon J5| this is a dc 
bias current). In that case, the soliton can be placed in 
an effective potential similar to that shown in Fig. 1. If 
it had lost all its kinetic energy [f22||6|, then we can use 
the static soliton as an initial condition. In fact, we think 
that the conditions for soliton tunneling could have been 
satisfied in the experimental situation described in [^2| . 

We conclude that the phenomenon of soliton tunneling 
is robust and generic. We believe this phenomenon can be 
observed also in other physical systems bearing solitons, 
topological defects, vortices, spiral waves, etc. 
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FIG. 1. Potential V(x) for the soliton escaping problem. 
The inset shows the force F(x). Note that in all figures the 
quantities plotted are dimensionless. 

FIG. 2. Numerical simulation of the soliton tunneling with 
sub-barrier kinetic energy. The pale curve is the potential and 
the bold curve is the soliton. The inflexion point is approx- 
imately the center of mass, (a) t = 0, V(t = 0) = 0; (b-f) 
show the dynamics in succesive time instants. 
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FIG. 3. Conditions for the existence of soliton tunneling 
with sub-barrier kinetic energy. The curve Fo = Fm 6 S sepa- 
rates the upper zone, where the soliton tunneling is possible 
from the lower zone, where the soliton tunneling is impossible. 
The filled circles represent numerical experiments. 

FIG. 4. Conditions for the existence of soliton tunneling 
with sub-barrier kinetic energy involving the parameters Fo 
and Vm- The other parameters remain fixed. 

FIG. 5. Disordered array of inhomogeneities. The soliton 
can move through the whole array. 

FIG. 6. Combination of impurities for Eq. (JsJ) . The soliton 
tunneling with sub-barrier kinetic energy is possible from zone 
B to zone C. 
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